IDENTIFICATION OF BOUNDARY CONDITIONS 
USING NATURAL FREQUENCIES 

AKHTYAMOV A. M."'i, MOUFTAKHOV A. V.^'^ 

° Department of Differential Equations, Bashkir State University, 
'\d '. Ufa, Russia. 

O ■ ^ Department of Mathematics and Statistics, Bar-Ilan University, 

^ ■ Ramat-Gan, Israel. 

M 

Abstract 

The present investigation concerns a disc of varying thickness of whose 
^r^ ■ flexural stiffness D varies with the radius r according to the law D = Dqt"^ , 

where Dq and m are constants. The problem of finding boundary conditions 

p, \ for fastening this disc, which are inaccessible to direct observation, from the 

C/^ ■ natural frequencies of its axisymmetric flexural oscillations is considered. The 

(-| . problem in question belongs to the class of inverse problems and is a completely 

"j^ \ natural problem of identification of boundary conditions. The search for the 

unknown conditions for fastening the disc is equivalent to finding the span of 

the vectors of unknown conditions coefficients. It is shown that this inverse 

problem is well posed. 
^ I Two theorems on the uniqueness and a theorem on stability of the solution 

^ ■ of this problem are proved, and a method for establishing the unknown condi- 

Q , tions for fastening the disc to the walls is indicated. An approximate formula 

t:J" \ for determining the unknown conditions is obtained using first three natural 

^ ' frequencies. The method of approximate calculation of unknown boundary 

Q . conditions is explained with the help of three examples of different cases for 

"^ \ the fastening the disc (rigid clamping, free support, elastic fixing). 

d • Keywords: Boundary conditions, a disc of varying thickness, inverse prob- 

lem, Plucker condition. 
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1. Introduction. Discs are parts of turbines and other various 
devices (see |I]-[6j). If it is impossible to observe the disc directly, the 
only source of information about possible defects of its fastening can be 
the natural frequencies of its flexural vibrations. The question arises 
whether one would be able to detect damage in disc fastening by the 
natural frequencies of its symmetric flexural vibrations. This paper gives 
and substantiates a positive answer to this question. 

The problem in question belongs to the class of inverse problems 
and is a completely natural problem of identification of the boundary 
conditions. 
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Closely related formulations of the problem were proposed in [3, [H]. 
Contrary to this, in this paper it is not the form of the domain or size 
of an object which are sought for but the nature of fastening. The 
problem of determining a boundary condition has been considered in 
[Hj. However, as data for finding the boundary conditions, we take not 
a set of natural frequencies, but not condensation and inversion (as in 

isi). 

Similarly formulated problems also occur in the spectral theory of 
differential operators, where it is required to determine the coefficients 
of a differential equation and the boundary conditions using a set of 
eigenvalues (for more details, see PIj-^^). However, as data for finding 
the boundary conditions, we take one spectrum but not several spectra 
or other additional spectral data (for example, the spectral function, 
the Weyl function or the so-called weighting numbers) that were used 
in PIIj-[IE]- Moreover, the principal aim there was to determine the 
coefficients in the equation and not in boundary conditions. The aim 
of this paper is to determine the boundary conditions of the eigenvalue 
problem from its spectrum in the case of a known differential equation. 

The problem of determining a boundary condition using a finite set 
of eigenvalues has been considered previously in |Tn|-[IB]- In contrast 
to these papers, we think it is necessary to determine not the type of 
fastening of a string, membrane or beam but the type of disc fastening 
of varying thickness. 

2. Formulation of the direct problem. The problem of axisym- 
metric flexural oscillations of disc of varying thickness without a hole in 
the centre is reduced to following differential equation ^ 

d'^w 1 dw\ dD f d'^w u dw\ 1 r ,d'^w 

r dr ) dr \ dr'^ r dr ) r dt"^ 




^1) 
where w is the deflection at radius r, i? is the flexural rigitiy, h is the 

thickness, p is the density and v is Poisson's ratio. 

\i D = Dor'^, where Dq = E hQ/12{l — z/^), and m are constants, 

the thickness h is given by h = ho r™/^. Then Equation ([T]) becomes 



r + r (2 m + 2) +r (m + z/m— 1 + m ^ 



w 



/f'Y*'^ //7^ f~i'Y'^ 



+r(l + z/m — lym — m) -;^ = ^ r 



dw ^ 12/j(l-Z/^) ^(i2_2m)/3 ^ 

6»r E/i2 5^2' 
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For vibrations, we write 



w = u(r)e'''^ 



and hence obtain the following equation for u{r): 

A (Jj il o /^ ^\ ^ tX rj ^ rj , (X LL 

r — — j- + r (z m + zj — -^ ^r [m + um — 1 + m j — -^ 
+r(l + V m^ — V m — m) — ^j.,'^'^ ^ ^(i2-2m)/3 ^ = q. 

Although it is impossible to solve the problem for arbitrary m exactly, 
some particular solutions in terms of Bessel functions are found for a 
number of positive values of m and corresponding values of Poisson 
ratio ly in [T^ . 

For the sake of being definite, we consider the particular case of 
Equation ^ with h = hor'^'^, Poisson ratio z/ = 1/9, 2f^jp — 1; ^^^ 
radius R = 1. 

The problem of axisymmetric flexural oscillations of a disc of varying 
thickness without a hole in the center for this particular case is reduced 
to following eigenvalue problem [H 



^d^u ^d^u 47 r'^ d'^ u Irdu ^ 8/3 

dr^ ^^3 9 (^^2 g ^^ 

Ui{u) = ai (Li u)^^^ + a4 {L4 u)^^^ = 0, (4) 

U2{u) = a2 {L2 u)^^^ + 03 (L3 u)^^^ = 0. (5) 

Here s^ = — ^'^Eh'^~^ ~ 21 eK^ ~ '^^ ^^ ^^^ eigenvalue parameter, E is 
the elasticity modulus, uo is the natural frequencies parameter, Ui{u), 
U2{u) are the linear forms which characterize conditions for fastening 
the plate to the walls (rigid clamping, free support, free edge, floating 
fixing, elastic fixing), 

, . duir) d^uir) 1 du 

Liu = u{r), L2U = —- — , Lsu = ——^ + -—--, 

dr dr^ \3r dr 



dr 



d'^u(r) 1 du 
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The solution to (El) has form [TH] 

•3^2/3^1/2- 



u = r 



-2/3 



Ci/i 



/3r2/3si/2x 

+ C2 Jl — ^ — + 



^c.Y,('-f^:^)+c,Yj'"-"^'"^ 



Here /i(-), Ji{-), Yi{-) are the conventional notations for first-order 
Bessel functions of real and imaginary argument [2D]. 

For a continuous plate (without a hole in the center), the constants 
C3 = C4 = 0. 

Instead of the (JH), (J5I), we can write the conditions 






i = 1, 2. 



(6) 



where an = ai, au = 0, aia = 0, au = 04, a2i = 0, 022 = 02, 023 = ^3, 
a24 = 0. 

We shall now formulate the direct eigenvalue problem (jHI, ©: it is 
required to find the unknown natural frequencies of the oscillations of 
the plate from the linear forms f/i(w), U2{u). 

The natural frequencies cjj are the corresponding positive eigenvalues 
of problem (EI), (EI) (see jHI). The non-zero eigenvalues of problem 
are the roots of the determinant 



Ms) 



Ui{ui) Ui{u2] 

U2{ui) U2{U2] 



= Mi2/l(s)+Mi3/2(s)+M24/3(s)+M34/4(s), 



where 

ui = wi(r, s) = r"^/^ /i (I r^/^ s^/^) , U2 = U2{r, s) = r"^/^ Ji (| r^/^ s^/^) 

are linearly independent solutions of Equation 



(7) 



/i(s) = [Liwi(r, s) • L2W2(r, s) - L2Wi(r, s) • Li W2(r, s) J^^^ , 
/2(s) = [Liwi(r, s) • L3W2(r, s) - L3Wi(r, s) • Li W2(r, s) J^^^ , 
fz{s) = [L^ui{r,s) ■ L2U2{r,s) - L2Ui{r,s) ■ L4 W2(r, s) J^^^ , 
/4(s) = [L4Wi(r, s) • L3W2(r, s) - L3Wi(r, s) • L4 W2(r, s) J^^^ . 

Thus, finding the three natural frequencies ui is equivalent to found- 
ing of three roots si of A(s). 





c = 10-5 


c= 1 


c= 10 


c= 102 


c = 10^00 




0.085457 
3.1122 
5.4634 


1.5178 
3.1145 
5.4651 


2.6517 
3.1561 
5.4813 


3.0663 
4.5575 
5.7412 


3.0739 
5.1995 
7.3054 



Table 1: Dependence c on Sj. 



For example, if boundary conditions (El) have the form 
Ui{u) = [cLiu{r) - L^u{r)l^, = 0, U2{u) = [L2u{r)l^, = 
(elastic clamping). For different c we have the different Sj. See Table 1. 

Thus, knowing c it is possible to find Si by standard methods. The 
solution to this direct problem presents no difficulties. The question 
arises whether one would be able to do the reverse and find c knowing 
Si. In a broader sense it may be stated as follows. Is it possible to derive 
unknown boundary conditions with a knowledge of s^? The answer to 
this question is given in the next section. 

3. Formulation of the inverse problem. The mathematical (di- 
rect) problem is an eigenvalue problem for a homogeneous linear fourth 
order equation, set up in the interval < r < 1, accompanied by two 
linear homogeneous boundary conditions at r = 1; the sought solutions 
must be bounded for r = 0. The boundary conditions depend on 4 
scalar coefficients (namely 3 because of the homogeneity). 

Now we shall formulate the inverse of eigenvalue problem (El), (El): 
it is required to find the unknown linear forms t/i(w), U2{u) from the 
natural frequencies of the disc oscillations. 

We shall denote the matrix, consisting of the coefficients aij of the 
forms Ui{u) and U2{u), by A and its minors by Mij-. 



A = 



ail CLi2 ai3 ai4 

a2l a22 CL23 0^24 



M 



«J 



an aij 
a2i a2j 



The search for the forms Ui{u) and U2{u) is equivalent to finding the 
span (ai, a2) of the vectors a^ = (an, 0^2, cia, ca^^ {i = 1, 2). 

Hence, in terms of eigenvalue problem (jSI), (El), the inverse problem 



constructed above should be formulated as follows: the coefficients aij of 
the forms Ui{u) and U2{u) of problem (EI), (E| are unknown, the rank of 



the matrix A to make up these coefficients is equal to two, the non-zero 
natural frequencies Sk of problem (El), (EI) are known and it is required 
to find the span (ai, a2) of the vectors a^ = {an, ai2, ai^, ai^)^ {i = 1,2). 

4. The uniqueness of the solution. Together with problem 
let us consider the following eigenvalue problem 



Ad'^u ^ od^u 47r^(i^M 7rdu o o/o 
"4 h 6 r \ s r ' u = 



+ 

4 






i = l,2, 
0, 622 



(9) 

&2, &23 = h, 



where 611 = 61, 612 = 0, 613 = 0, bu = &4, &21 
624 = 0. 

We denote the matrix composed of the coefficients bij of the forms 
Ui{u) and U2{u) by B and its minors by Mij-. 



B = 



bn 612 &13 bu 

621 &22 &23 ^24 



M,j = 



bu bij 
b2i b2j 



The span of the vectors b^ = (6^1, 6^2, &i3, ^^4)^ (^ = 1, 2) is denoted 
by (bi,b2). 

Theorem 1 (on the uniqueness of the solution of the inverse problem) . 
Suppose the following conditions are satisfied 



rank74 = rank 5 = 2. 



(10) 



If the non-zero eigenvalues [sk] of problem (EI), (EI) and the non-zero 
eigenvalues {sk] of problem (El), (ED are identical, with account taken 
for their multiphcities, the spans (ai, a2) and (bi, b2) are also identical. 
Proof. The non-zero eigenvalues of problem (EI), (ED are the roots of 
the determinant 



Ms) = 



Ui{ui) U\{U2) 

U2{ui) U2{U2) 



= Mi2/i(s)+Mi3/2(s)+M24/3(s)+M34/4(s). 



In additional to the roots identical to the non-zero eigenvalues of the 
problems, the determinants A(s) and A(s) can also have the root s = 
of finite multiplicity. 



Since A(s) ^ 0, A(s) ^ are entire functions in s of order 1/2, it 
follows from Hadamard's factorization theorem (see [21]) that determi- 
nants A(s) and A(s) are connected by the relation 

A{s) = Cs''A{s), 

where A; is a certain non-negative integer and C is a certain non-zero 
constant. From this we obtain the identity 

(Mi2 -Cs'^Mu) /i(s) + (Mi3 - Cs^ M13) /2(s)+ 

+ (M24-CS^M24)/3(S) + (M34-CS^M34)/4(S)=0. ^ ^ 

Note that the number k in this identity is equal to zero. Actually, let 
us assume the opposite: k ^0. Then the functions fi{s) {i = 1,2,3,4) 
and also the same functions multiplied by s^ are linearly independent. 
Indeed, using MAPLE, we get 

27 2 729 4 6561 ^ 177147 o ^, ,„, 

f. s) = s -\ s s -\ s + Ofs 

■^ ^ ^ 64 131072 671088640 38482906972160 ^ ' 

^ , , 3 2 567 4 9477 ^ 373977 o ^, .n, 

■^ ^ ^ 16 32768 167772160 9620726743040 ^ ^' 
. / N 81 4 5103 . 216513 „ ^, ^n, 

•^^^ ^ 1024 5242880 150323855360 ^ ^' 

/4(S) = A54_^^^6^^5309_^3 10)_ 

•^ ^ ^ 128 32768 2684354560 ^ ^ 

The determinant of the matrix 4 x 4 of the coefficients at s* {i = 
2, 4, 6, 8) in these representations of the functions fj [j = 1,2,3,4) is 
not equal 0. It now follows that eight functions fj, s'^ fj (j = 1,2, ... , 4) 
are linearly independent for all k > 8. 

Further, the determinant of the matrix 8 x 8 of the coefficients at 



s* (i = 4, 6, . . . , 18) in the series expansion of the functions fj, s f. 
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(j = 1, 2, . . . , 4) is not equal under k— 2, k— 4 or k=6. This means that 
eight functions fj{s), s^ fj{s) {j = 1,2,..., 4) are hnearly independent 
under k=2, k=4 or k=6. 

This completes the the proof of linear independence of the functions 
/„s^/, (j = l,2,...,4)foralU7^0. 

From this and identity f pTl ) we obtain 



M12 = Mi3 = M24 = M34 = M12 = Mi3 = M24 = M34 = 0, 



which, in combination with Mu = M23 = 0, contradict condition I HM 
of the theorem. 

Hence, k = 0. From this and from identity ( TTTl l. by virtue of the 
hnearly independence of the corresponding functions, we obtain 



(M12, Mi3, Mi4, M23, M24, MuY = C (M12, Mi3, Mi4, M23, M24, M34)' , 

which is equivalent to the proportionahty of the bivectors ai A a2 and 
bi A b2. 

It is well-known ^j that there is a natural one-to-one correspondence 
between the classes of non-zero, proportional bivectors and the two- 
dimensional subspaces of a vector space. In this correspondence, a vector 
product X1AX2 of the vectors of its arbitrary basis xi and X2 corresponds 
to each subspace and a subspace (xi,X2) corresponds to each bivector 
Xi A X2. It therefore follows from the last equation that (ai, a2) which 
it was required to prove. 

5. Exact solution. It has been shown above that the problem 
of finding the unknown linear forms Ui{u) and U2{u) from the natu- 
ral frequencies of axisymmetric flexural oscillations of a disc of varying 
thickness has a unique solution (in the sense that the spans, composed 
of the coefficients of these linear forms, are uniquely defined). The next 
question is how this solution can be constructed. 

This section deals with solving this problem and constructing exact 
solution by the first three natural frequencies cui. 

Suppose si, S2, S3 are the values corresponding to the first three 
natural frequencies cUi. We substitute the values si, S2, S3 into A(s) 
and obtain a system of three homogeneous algebraic equations in the 
four unknows M12, M13, M24, M34 

Mu /i(s,) + Mi3 /2(s,) + M24 /4(s,) + M34 /4(s,) = 0, ^ = 1, 2, 3. 

(12) 
The resulting set of equations has an infinite number of solutions. 
If the resulting set has a rank of 3, the unknown minors can be found 
in accurate to a coefficient. The unknown span and its basis are found 
from the minors using well-known methods of algebraic geometry [23j. 
We also note [12] and f2^, in which detailed method for solution of this 
problem was given. 



The resulting set of equations has an infinite number of solutions. It 
follows from the uniqueness theorem already been proved that the un- 
known minors M12, M13, M24, M34 can be found accurate to a constant. 
Hence, the resulting system must have a rank of 3 and the solution is 
determined accurate to a constant multiplier. 

If the minors M12, M13, M24, M34 are found accurate to a constant 
then the unknown span and its basis are found from the minors using 
well-known methods of algebraic geometry [23j. We also note [12] and 
[2l] . in which a detailed method for solution of this problem was given. 

For example, if M12 = 1, then 



A = 



10 -M24 
1 Mi3 



if Mi3 = 1, then 



A = 



1 






M12 





1 



-M34 





(13) 



(14) 



These reasons prove 

Theorem 2 (on the uniqueness of the solution of the inverse problem). 
If the matrix of system ( TT^ has a rank of 3, the solution of the inverse 
problem of the reconstruction boundary conditions (HI), (0) is unique. 

Remark. Theorem 2 is stronger than theorem 1. Theorem 2 use 
only three natural frequencies for the reconstruction of boundary con- 
ditions and not all natural frequencies as in theorem 1. It follows from 
the uniqueness theorem 1 already proved that the unknown minors 
^12,^13,^24,^34 can be found accurate to a constant by all natu- 
ral frequencies. To prove that the unknown minors M12, M13, M24, M34 
can be found accurate to a constant by three natural frequencies is very 
easy problem in concrete cases. But in the common case to prove that 
set in^ has a rank of 3 is very failed to do this. Yet our uniqueness 
theorem 1 suggests this result in the common case. 

6. Approximate solution. Since small errors are possible when 
measuring natural frequencies, the problem arises of finding an algo- 
rithm for the approximate determination of the type of fixing disc from 
the first three natural frequencies found with a certain error. 

If the values Sj {i = 1, 2, 3) are approximately the same as the first 
three exact eigenvalues of problem (EI), (EI), and set ( 0^ has a rank 



equal to three, the unknown minors can be determined accurate to a co- 
efficient. Further, the problem arises of finding the unknown span from 
the values M12, M13, M24, M34. However, the values M12, M13, M24, M34 
cannot be minors of a matrix. So this problem is not trivial. We must 
find minors P12, P13, P245 ^34 close to the values M12, M13, M24, M34. 
This problem is solved with the help Lagrangian multiplier method and 
algebraic geometry. 

It is known from algebraic geometry that the numbers P12, P13, Pu, 
P23, P24, Pm are minors of some matrix, if and only if the following 
condition is satisfied 

P12P34 - P13P24 + P14P23 = 0. 

This condition is called Plucker condition (see [^1^ ). 
Using Pu = P23 = 0, we get 

P12P34 - PuP24 = 0. (15) 

By definition, put xi = P12, X2 = P34, x^ = P13, 0:4 = — P24- Using 
this definition, we get Plucker condition 

xi 2:2 + X3 2:4 = 0. (16) 

It characterizes a surface in the 4-dimensional space. 

By definition, put yi = Mu, 2/2 = M34, 7/3 = M13, 7/4 = -M24. 

If Plucker condition for these numbers is realized, then M12, M34, 
Mi3, M24 are minors of some matrix and corresponding boundary con- 
ditions are found from the minors using the methods in section 4. 

If Plucker condition for numbers M12, M34, M13, M24 is not realized, 
the required minors P12, P13, P24, P34 are found with the help Lagrangian 
multiplier method and algebraic geometry. 

Indeed suppose P(a:i,a:2, X3, X4,j9) is Lagrange function 

{xi - yif + {x2 - 2/2)^ + {xi - ysf + (2:4 - y^f + 2;? (xi X2 + X3 X4), 

where 2j9 is the Lagrangian multiplier. If we find minimum of the 
function F{xi, X2: X3, X4, jy), then we obtain the minors P12, P13, P24, P34 
most close to the values M12, M13, M24, M34. 
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The minimum of the function F(xi, 0:2, 0:3, 0:4, ]?) is found from equa- 
tions 



(17) 



F'^^{xi,X2,Xi,XA,p) = 2(a:i - yi) + 2px2 = 0, 

F^^{xi,X2,X3,X4,p) = 2(2:2-2/2) + 2^3X1 = 0, 
F^^{xi,X2,X3,X4,p) = 2(2:3 -7/3) + 2^3X4 = 0, 
F^^(xi,X2,X3,a:4,;p) = 2(x4-?/4) + 2;pX3 = 0, 

F2p(a:i,a:2,a:3,a:4,j3) = xi 0:2 + 2:3 0:4 = 0. (18 



By definition, put 

X = (xi, X2, X3, X4), y = {yu 2/2, ?/3, ?/4), 

X* = (X2, Xi, 3:4, X3), Y* = (?/2, 2/1, 2/4, 2/3), 

(X, >^) = 3:1 2/1 + 0:2 2/2 + 3:3 2/3 + 3:4 2/4- 

With the preceding notations Equations f fTTl l. f p^ are identical to 
the following equations 

Y = X+pX\ (19) 

(X,X*) = 0. (20) 

It follows from f HUl l and (EOl) that the vector X* is orthogonal to the 
vector X, and the vector X = (xi, 2:2, X3, X4) is orthogonal projection 
of the vector Y = (2/1, 2/2, 2/3 5 2/4) on the surface ( TT^ . 

Having solved f H^ as set of linear equations with the unknowns 

xi, 0:2, 3:3, X4, we get 

-^ = Y^(>'-pn. (21) 

From i HTh it is easy to get 

X' = -^{Y'-pY). (22) 

1 — J9^ 

Substituting (EII) for X and (El) for X* in (EDI), we obtain 

{Y -pY', Y* -pY) = 0. 
Notice that 

(y, y*) 7^ 0, (y*, y*) = (y, y), (y*, y) = (y, y*). 



11 



Therefore, 



,2 o „ (5^. y) 



^-2^(Tfy^ + ^ = °- 



This quadric equation has two roots 

{Y,Y)TViY^Yy-(y^y*y 



If X is close to Y, then |j3| << 1 and thus we have 

P- ^Y,Y*) • ^^^^ 

The vector X can be found using (J2II) and (J23I)- The coordinates 
Pi2: Pu: P2A1 Pm of X are minors of a matrix. This matrix and corre- 
sponding boundary conditions are found from the minors P12, P13, P24, ^34 
using the methods in section 4. 

7. Stability of the solution. In this section we study continuity 
of the solution of the inverse problem with respect to Sj. It is shown 
that small perturbations of eigenvalues Sj [i = 1,2,3) lead to small 
perturbations of the boundary conditions. 

Let Si [i = 1,2,3) be eigenvalues of problem (El,®, and s^ {i = 
1, 2, 3) are such values that 

\si — Si\ < 6 << 1, i = 1,2,3; 
and R is such a number that 

\si\ < R, |si| < R- 

Theorem 3. Suppose that one of the third-order minors of the ma- 
trix ||/fc(si)||3x4 is substantially non-zero. If |si — s^l < (5 << 1, then 
the boundary conditions of problem (EI),© are close to the boundary 
conditions of problem (ED,(|E1)- 

Proof. By virtue of continuity of the functions (C|), we have 

\fk(S^) - fM)\ < M 5, ^ = 1,2,3, /c = 1,2,3,4, (24) 
where M = max M^, Mh = max fkiz). 

fc=l,2,3,4 \z\<R 
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First, let us prove that 

\M,j - M,j\ < Ci6. 

Recall that one of the third-order minors of the matrix ||//c(si)||3x4 is 
substantially non-zero by the hypotheses of Theorem 3, and so 

rank||/fc(s^)||3x4 = 3, z = l,2,3. A; =1,2, 3, 4. (25) 

To be precise, assume that 

mod (det||/fc(s,)||3x3) »0, z= 1,2,3, /c = 1,2,3. (26) 

It follows from (El) that 

mod (det||/fc(?,)||3x3) »0, z = l,2,3, A; = 1,2,3. (27) 

It can be shown by direct calculations that set dT^ has the solution 



:28i 



Like previously, the set 

Mu /i(s.) + Mis f2{s^) + M24 /3(s,) + M34 h{s,) = 0, z = 1, 2, 3. 

(29) 
has solution 



Mi2 = det( fk{si) 3x3), 


z = l,2,3. 


A; = 2, 3, 4, 


Mi3 = det( fk{s^) 3x3), 


z = l,2,3. 


k = 1,3,4, 


M24 = det( fk{s^) 3x3), 


z = l,2,3. 


A: = l,2,4, 


M34 = det( fk{si) 3x3), 


z = l,2,3. 


A; = l,2,3. 



Mi2 = det( fk{s,) 3x3), 


z = l,2,3. 


A; = 2, 3, 4, 


Mi3 = det( fk{s^) 3x3), 


z = l,2,3. 


A; = l,3,4, 


M24 = det( fk{s,) 3x3), 


z = l,2,3. 


A:=l,2,4, 


M34 = det( fkis^) 3x3), 


z = l,2,3. 


k = 1,2,3. 



(30) 



By definition, put 

\k = fk{Si) — fk{Si), Afc = ||Ajfc||3xi, Fk = \\fk{Si)\\3xl, 

{i = 1,2,3, k = 1,2,3,4). 
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In the new notations, we have 



|Mi2 - M12I = mod (det(||/fc(si)||3x3) - det(||/A:(si)||3x3)) = 

= mod (det||Fi + Ai, F2 + A2, F3 + A3II - det ||Fi, F2, F3II) = 
= mod (det||Fi, F2 + A2, F3 + A3II + 

+ det||Ai, F2 + A2, F3 + A3II - 

-det||Fi, F2, F3II) = 
= mod (det||Fi, F2, F3II + det ||Ai, F2, F3II + 

+ • • • + det II Ai, A2, A3II - det ||Fi, F2, F3II) = 
= mod (detlJAi, F2, F3II + • • • + det ||Ai, A2, A3II). 

It follows from (El) that 

\^^k\<M5. (31) 

Calculating the determinant and using (Ell), we get 

mod (det ||Ai, F2, F3II) < 6 • M^ • (5. 

Similarly, mod (det ||Ai, A2, F3II) < 6 • M^ • (5^ < 6 • M^ • (5. 
Arguing as above, we see that 

mod (det ||Ai, F2, F3II + h det ||Ai, A2, A3II) < 7 • 6 • M^ • (5. 

By definition, put Ci = 7 • 6 • M^ . Then, 

\Mi2-Mi2\= mod (det(||/fc(i'OI|3x3)-det(||/fc(sO||3x3)) <Ci5. 
Continuing this line of reasoning, we see that 

\M,j - M,j\ < Ci5. 
This implies that 



\Y -Y\ = 



\ "^{yk-Vky <'2Cid, 
\ fc=i 



(32) 



where 



Y = (Mi3, -M24, Mi4, M 



23^ 



\T 



Y = (Mi3, -M24, Mi4, M23) 



T 



14 



By O, (EZD this means that Y is clos^to Y. _ 

Now let us prove that X = (P12, Pm, Pi3, —P24)^ close to 

X = (^12, ^34, ^13) —^24) • 

Combining (J23I) and (E2D, we see that 



P 



{Y ,Y) - ^ {Y ,YY - {Y ,Y*Y 



{Y,Y^ 



close to p. This means that X = (F - pY*)/{l - p) close to (E 
Finally, let us prove that the boundary conditions of problem (jHI) 
are close to the boundary conditions of problem (EI),®, where 

\si — Si\ < 5 « 1. 

We know already (see [21]) that the coefficients of the boundary 
conditions are the minors of matrix A as in JH^ . This means that if X 
is close to X, the boundary conditions of problem (jHl),® are close to 
the boundary conditions of problem (EI),®- 

This completes the proof. 

Computer calculations confirm the stability of the solution of the 
inverse problem. The order of error often hardly different from the error 
in the closeness of values of Aj and A^ and only in some cases it can be 
deteriorated by four orders of magnitude. So the measurement accuracy 
of instruments to measure natural frequencies must exceed accuracy to 
measure boundary conditions by four orders of magnitude. 

It follows from theorem 2 and 3 that the inverse problem is well 
posed, since its solution exists, is unique and continuous with respect to 
Sj {i = 1,2,3). 

8. Examples. We use dimensionless variables in the numerical 
examples. 

Example 1. If si^''^ = 3.0739, S2^''^ = 5.1995, F3 ^^^ = 7.3054 
correspond to the first three natural frequencies Ui determined using 
instruments for measuring the natural frequencies with an accuracy of 
10^^, then the solution of set dT^ . accurate to a constant, has the form 

M12 = 645330 C, Mi3 = 19.947 C, M24 = 2.4157 C, M34 = C. 
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Pi3 = 19.94703753 C, 
P34 = 0.000749141 C. 



Using (EH) and (ESJ, we get 

P12 = 645330 C, 
P24 = 2.416009841 C, 

Suppose C = I/P12; then from f H^ . we obtain 

1 -3.75565 • 10"^ 
1 0.00031 



A 



1/2 _ 



1/2 _ 



1/2 _ 



Note that the numbers si ' = 3.0739, S2 = 5.1995, S3 
7.3054 presented above are almost the same as the first three exact 
values si^/^, 52^^^, 53^'^ corresponding to rigid clamping. This means 
that the unknown disc fastening inaccessible to direct observation has 
been correctly determined. 

Example 2. If si ^/^ = 1.8312, S2 ^'^ = 4.4629, S3 ^/^ = 6.6502 
correspond to the first three natural frequencies coi determined by a 
frequency meter with an accuracy of 10~^, then the solution of set ( TT^ . 
accurate to a constant, has the form 

M12 = -9.31 C, Mi3 = 201420 C, M24 = -2.6702 C, M34 = C. 

Using (EH and O, we get 

P12 = -9.310013258 C, P13 = 201420 C, 
P24 = -4.62276- 10-^ C, P34 = 1.00012342 C. 

Suppose C = I/P13; then from I HM . we obtain 

1 -4.96536 • 10"^ 



A = 



-4.62219-10-^ 1 



1/2 ~ 1/2 







1/2 



Note that the numbers si , S2 , S3 ' presented above are almost 
the same as the first three exact values si^'^, S2^'^, S3^/^ corresponding 
to the free support. This means that the unknown disc fastening inac- 
cessible to direct observation has been correctly determined. 

Example 3. If si^''^ = 1.5178, S2 ^^^ = 3.1145, S3 ^^^ = 5.4651 
correspond to the first three natural frequencies ui determined by means 
of instruments for measuring natural frequencies with an accuracy of 
10^^, then the solution of set (O, accurate to a constant, has the form 



Mi2 = -140260C, Mi3 = 154.74C, M24 = 140150C, M34 = C. 
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Using (EH) and (ESD, we get 



Pi2 = -140260 C, 
P24 = 140150 C, 



p^3 = 76.931 C, 
P34 = -76.870 C. 



Suppose C = I/P12; then from ( TTTIl ). we obtain 



1 
1 



-0.99922 

0.00054849 



Note that the numbers si , S2 , S3 presented above are al- 
most the same as the first three exact values si"^'^, S2^'^, 53^/^, which 
correspond to elastic fixing with matrix 



A = 



10 0-1 
10 



This means that the unknown disc fastening inaccessible to direct ob- 
servation has been satisfactorily determined. 

Thus, the form of the disc fastening of varying thickness can be de- 
termined from the first natural frequencies measured by special instru- 
ments. 

Note that we choose such a particular variation law for the thickness 
of the disc to be precise. If the variation law for the thickness of the 
disc is different from that adopted in this study, then the mathematical 
formulation and the proposed procedures of solving of the direct and 
inverse problems remain valid. In this case we must substitute linearly 
independent solutions of corresponding differential equation for ui and 
W2 in (0), whose are not singular at r = 0, and the corresponding linearly 
forms for Li, L2, L3, L4 in (El). 

Direct problems on hydroelasticity and aeroelasticity are considered 
in PSl [2Ej- Similar inverse problems can be solved by means of the 
method proposed in this paper. 
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